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Abstract. We report the current status of our systematic studies of the QCD 
thermodynamics by lattice QCD simulations with two flavors of improved Wilson 
quarks. We evaluate the critical temperature of two flavor QCD in the chiral 
limit at zero chemical potential and show the preliminary result. Also we 
discuss fluctuations at none-zero temperature and density by calculating the quark 
number and isospin susceptibilities and their derivatives with respect to chemical 
potential. 



1. Introduction 

In order to extract unambiguous signals for the QCD phase transition from the 
heavy-ion collision experiments, quantitative calculations from first principle are 
indispensable. At present, the lattice QCD simulation is the only systematic method 
to do so and various interesting results have been already reported. So far, most 
lattice QCD studies at finite temperature (T) and chemical potential (/i q ) have been 
performed using staggered-type quark actions with the fourth-root trick of the quark 
determinant. Thus, the other actions such as the Wilson-type quark actions are 
necessary to control and estimate systematic errors due to lattice discretization. 

Several years ago, the CP-PACS Collaboration has studied the QCD 
thermodynamics using the Iwasaki (RG) improved gauge action and the two-flavor 
clover improved Wilson quark action pQ. We revisit this action armed with recent 
techniques at finite chemical potential [i q . In this report, we present on our preliminary 
results for the critical temperature at \i q = as well as the quark number and isospin 
susceptibilities at finite fi q . The latter are related to the physics of the possible critical 
point in the (T, plane. 
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2. Critical temperature of two-flavor QCD 

The critical temperature (T c ) is one of the most fundamental quantities in the QCD 
thermodynamics. We renew the analysis of T c done in [I] performing additional 
simulations at p = 6/g 2 = 1.8 for Nt = 4 and 1.95 for N t — 6. We determine 
the pseudo-critical point P pc defined from the peak of Polyakov loop susceptibility on 
V 3 x N t — 16 3 x 4 and 16 3 x 6 lattices, as a function of the hopping parameter K. 

It is confirmed in Ref. pQ that a subtracted chiral condensate satisfies the scaling 
behavior with the critical exponents and scaling function of the three-dimensional 0(4) 
spin model. The reduced temperature t and external magnetic field h are identified 
to t ~ P — Pd and h ~ m q , respectively, where Pd is the critical transition point in 
the chiral limit. Assuming that the pseudo-critical temperature from the Polyakov 
loop susceptibility follows the same scaling law as the 0(4) spin model, t pc ~ h y with 
y = 0.537(7), we fit the data of P pc by P pc = Pet + Ah 1 ^ with two free parameters, p c t 
and A, in the range of P = 1.8-1.95 for N t = 4 and P = 1.95-2.10 for N t = 6. First, 
we adopt the definition m q a ~ 1/K — 1/K C as the quark mass where K c is the chiral 
point where the pion mass vanishes at T = for each p. See Fig. Q] (left). We calculate 
the critical temperature T c in the chiral limit using T — l/N t a. The lattice spacing 
a is estimated from the vector meson mass assuming my(T = 0) = m p = 770 MeV 
at P c t on K c . By this procedure, we obtain the preliminary results T c = 183(3) MeV 
for Nt = 4 and 174(5) MeV for N t = 6. We also calculate p c t using the relation of 
m^ WI oc rripg, where mps is the pseudo-scalar meson mass and m^ WI is the quark mass 
obtained from the axial vector Ward-Takahashi identity. The results of T c are 173(3) 
MeV (N t = 4), 167(3) MeV (N t = 6) for h = (m PS a) 2 and 176(3) MeV (Nt = 4) 
for h = m^ WI a. We note that these 0(4) fits reproduce the data of P pc much better 
than a linear fit P pc = p ct + Ah. From these analyses, we tentatively conclude that 
the critical temperature in the chiral limit is in the range 170-186 MeV for N t = 4 
and 164-179 MeV for Nt = 6. There is still a large uncertainty from the choice of 
the fit ansatz. To remove this, we are performing further simulations at lighter quark 
masses. 

Next, we compare our results with those of the staggered quark action. We 
plot the results of the pseudo-critical temperature (T pc ) in unit of Sommer scale (tq) 
as a function of mps^o m Fig. Q] (right) together with those by the RBC-Bielefeld 
Collaboration using 2+1 flavor p4-improved staggered quark action [2]. As seen in 
this figure, results of T pc obtained by different quark actions seem to approach the 
same function of mpsro as Nt increases. 

A care is in order when we convert T c r$ (T pc ro) to T c (T pc ) in MeV using a 
physical value of tq . Because the phenomenological estimate of ro has large theoretical 
uncertainties, it looks convenient to adopt a lattice result. Unfortunately, lattice 
results for ro suffer from sizable ambiguities yet. The RBC-Biclcfcld Collaboration 
adopted the value ro = 0.469(7) fm, which was obtained by Gray et al. [3] from a 
bottomonium mass splitting using the AsqTad-improved staggered quark action. On 
the other hand, the CP-PACS+JLQCD Collaboration found r = 0.516(21) fm from 
a study of light hadron spectrum using the clover-improved Wilson quark action [4] . 
This leads to about 10 % difference in the value of T c and makes it difficult to naively 
compare results from different groups. 
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Figure 1. Left: The pseudo-critical point f3 pc as a function of m q ~ 1/K — 1/K C 
for Nt = 4 (circle) and Nt = 6 (square). Right: Comparison of T pc scaled by ro 
between the staggered quark action (open symbol) and the Wilson quark action 
(filled symbol) for Nt = 4 and 6. 

3. Hadronic fluctuations at finite /i q 

Hadronic fluctuations at finite density are observables closely related to the critical 
point in the (T,fj, q ) plane and may be experimentally detected by an event-by-event 
analysis of heavy ion collisions. The fluctuations can also be studied by numerical 
simulations of lattice QCD calculating the quark number and isospin susceptibilities, 
X q and xi- They correspond to the second derivatives of the pressure with respect 
to \i q and hi, where \ii is the isospin chemical potential. From a phenomenological 
argument in the sigma model, Xq is singular at the critical point, whereas xi shows 
no singularity there. 

We perform simulations at mps/mv = 0-65 and 0.80 on a 16 3 x 4 lattice. We 
calculate the Xq and Xi and their second derivatives with respect to fiq and fii at 
jtx ? = Hz = 0. (Note that the odd derivatives are zero at \i q = 0.) The details of the 
calculations are reported in Ref. [5] . 

The left panel of Fig.[5]shows Xq/T 2 (circle) and Xi/T 2 (square) at mps/my = 0.8 
and /i ? = /i/ = as functions of T/T pc . We find that Xq/T 2 and Xi/T 2 increase 
sharply at T pc , in accordance with the expectation that the fluctuations in the QGP 
phase are much larger than those in the hadron phase. Their second derivatives 
d 2 (Xg/T 2 )/d(p q /T) 2 and d 2 { X i /T 2 )/d(n q /T) 2 are shown in Fig. H (right). We find 
that basic features are quite similar to those found previously with the p4-improved 
staggered fermions [6]. d 2 (xi /T 2 )/d(fi q /T) 2 remains small around T pc , suggesting 
that there are no singularities in xi at non-zero density. On the other hand, we expect 
a large enhancement in the quark number fluctuations near T pc as approaching the 
critical point in the (T, fi q ) plane. The dashed line in Fig. [2] (right) is a prediction from 
the hadron resonance gas model, d 2 Xq/d/i 2 ~ 9x q /T 2 . Although current statistical 
errors in Fig. [2] (right) are still large, we find that d 2 (xq/T 2 )/d(/j, q /T) 2 near T pc is 
much larger than that at high temperatures. At the right end of the figure, values 
of free quark-gluon gas (Stefan-Boltzmann gas) for AT t = 4 and for N t = oo limit are 
shown. Since the lattice discretization error in the equation of state is known to be 
large at Nt = 4 with our quark action, we need to extend our study to larger Nt for 
the continuum extrapolation. 
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Figure 2. Left: Quark number (circle) and isospin (square) susceptibilities at 
= 1^1 = 0- Right: The second derivatives of these susceptibilities. 

4. Conclusion 

We reported the current status of our study of QCD thermodynamics using the two- 
flavor improved Wilson quark action. The critical temperature is estimated from the 
Polyakov loop susceptibilities on 16 3 x4 and 16 3 x 6 lattices. We discussed uncertainties 
from the chiral extrapolation and the scale parameter. Our preliminary result of T c 
in the chiral limit is in the range of T c = 170-186 MeV for N t = 4 and T c = 164-179 
MeV for N t = 6. 

The fluctuations of quark number and isospin densities were also discussed. 
Although the statistical errors are still large, we find that Xq seems to increase rapidly 
near T pc as [i q increases, whereas the increase of xi is n °t large near T pc . These 
behaviors qualitatively agree with the previous results obtained with the p4-improved 
staggered fermions. 

We also studied the static quark free energies, Debye screening masses and spatial 
string tension at finite temperature with the same action, which are reported in detail 
in Ref. [7]. 
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